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KINETIC  STUDY  OF  THE  COMBUSTION  OF 
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The  thermochemistry  and  kinetics  of  organophosphorus  compounds  were  studied,  with  BAC-MP4 
method  estimations  as  a  basis  and  PM3  semi-empirical  estimations  for  many  new  compounds.  New  group 
additivity  values  were  proposed  for  enthalpies  of  formation  at  298  K,  and  entropies  and  heat  capacities  of 
species  involving  pentavalent  phosphorus  bonded  to  carbon,  hydrogen,  oxygen,  fluorine,  nitrogen,  and 
sulfur  atoms.  The  kinetic  features  of  unimolecular  elimination  were  investigated  by  modeling  pyrolysis 
experiments  of  dimethyl  ethyl  phosphonate  (DEMP),  triethyl  phosphate  (TEP),  and  di-isopropyl  methyl 
phosphonate  (DIMP).  Rate  constants  were  proposed  for  four-  and  six-center  eliminations.  A  new  com¬ 
bustion  mechanism  was  developed  in  a  systematic  and  comprehensive  way  and  involved  41  phosphorus 
compounds  in  202  reactions.  All  possible  intermediates  were  taken  into  account.  This  reaction  mechanism 
was  applied  to  the  modeling  of  H2/02  flames  doped  with  dimethyl  methyl  phosphonate  (DMMP)  and 
trimethyl  phosphate  (TMP),  which  are  chemical  warfare  surrogates.  The  promoting  effect  of  the  agent 
and  the  species  profiles  were  correctly  reproduced.  The  decomposition  reaction  channels  were  studied; 
they  included  radical  reactions  and  molecular  eliminations.  The  cause  of  the  increase  of  the  reactivity  of 
the  flame  was  attributed  to  some  radical  combination  cycles  leading  to  an  increase  of  heat  release,  which 
enhances  the  global  reaction  rate  in  spite  of  the  radical  concentration  drop. 


Introduction 

Combustion  interest  in  thermochemistry  and  re¬ 
action  mechanisms  of  organophosphorus  com¬ 
pounds  derives  from  their  roles  in  incineration  of 
pesticides,  as  chemical  warfare  agents,  in  fire  sup¬ 
pression  [1-4],  and  in  potential  catalytic  applications 
in  aircraft  turbines  [5,6], 

Detailed  chemical  mechanisms  must  be  devel¬ 
oped  for  these  applications.  Few  experimental  stud¬ 
ies  have  been  reported.  In  pyrolysis,  below  1000  K, 
L’Homme  et  al.  [7]  determined  products  of  the  re¬ 
action  of  alkyl  phosphates;  Durig  et  al.  [8]  compared 
reactivity  of  six  compounds  and  proposed  a  decom¬ 
position  mechanism  by  unimolecular  eliminations. 
Zegers  and  Fisher  studied  thermal  decomposition  of 
diethyl  methyl  phosphonate  (DEMP)  [9],  triethyl- 
phosphonate  (TEP)  [10],  and  diisopropyl  methyl 
phosphonate  (DIMP)  [11]  in  the  range  700-900  K. 
They  measured  the  rate  of  unimolecular  reaction  of 
reactants  and  proposed  mechanisms  based  on  suc¬ 
cessive  molecular  eliminations.  In  combustion  con¬ 
ditions,  Werner  and  Cool  [12]  and  Korobeinichev  et 
al.  [13]  studied  the  addition  of  dimethyl  methyl 
phosphonate  (DMMP)  in  a  H2/02  flame;  both  mea¬ 
sured  phosphorus-containing  intermediates  and 
products.  Addition  of  trimethylphosphate  (TMP)  to 
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a  H2/02  flame  [14,15]  and  a  CH4/02  flame  [4]  was 
also  reported  by  Korobeinichev  et  al.,  with  a  quan¬ 
titative  determination  of  the  final  phosphorus  spe¬ 
cies  [4,15,16], 

Twarowski  proposed  a  detailed  mechanism  to  ex¬ 
plain  phosphine  addition  on  the  rate  of  recombina¬ 
tion  of  H  and  OH  [17],  including  6  species  and  33 
reactions.  The  non-tabulated  thermochemical  data 
were  evaluated  by  ab  initio  calculations.  An  ex¬ 
tended  version  [18]  was  proposed  involving  17  spe¬ 
cies  and  162  reactions,  reduced  after  analysis  to  35 
reactions.  Werner  and  Cool  proposed  a  19-step 
mechanism  for  combustion  of  DMMP  [12],  using 
Twarowski’s  data  for  the  small  species.  Korobeini¬ 
chev  et  al.  [4]  used  a  mechanism  with  34  reactions 
to  describe  TMP  destruction  in  flames  by  radical  re¬ 
actions  without  molecular  elimination. 

The  present  work  reports  evaluations  of  thermo¬ 
chemical  data  as  group  contributions,  allowing  esti¬ 
mation  of  thermochemical  properties  for  numerous 
new  species.  A  kinetic  study  of  the  decomposition 
mechanism  of  alkyl  phosphate  and  phosphonate  was 
performed  in  pyrolysis  and  in  flame  to  model  the 
effect  of  organophosphorus  additives  in  combustion 
and  the  chemistry  controlling  incineration  of  chem¬ 
ical  warfare  surrogates. 

Thermochemistry 

The  thermochemistry  of  phosphorus  compounds 
is  mostly  unknown.  Without  experimental  kinetic 
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TABLE  1 

Standard  enthalpy  of  formation  of  selected  compounds  (kcal  mol-1) 


Species 

This  Study 

BAC-MP4  [21] 

PM3  Calc. 

Benson  [19] 

PO(CH3)(OCH3)2 

- 196.4 

- 196.4 

- 189.2 

-211.1 

PO(H)(CH3)(OCH3) 

-131.1 

-131.5 

-134.1 

- 

PO(OH)(CH3)(OCH3) 

-  204.4 

-206.7 

- 195.6 

-225.4 

PO(OH)(OCH3)2 

-255.1 

-255.6 

-  250.6 

-271.2 

PO(OH)3 

-271.1 

-272.3 

-  270.2 

-299.6 

PO(OH)2(CH3) 

-212.4 

-210.9 

-209.1 

-239.6 

data  available,  a  postulated  kinetic  model  is  strongly 
dependent  on  estimation  of  the  heats  of  formation 
and  entropies  of  the  species,  Benson  [19]  proposed 
group  additivity  values  for  estimation  of  enthalpies, 
but  not  enough  to  allow  characterization  of  all  spe¬ 
cies  involved  in  a  detailed  mechanism.  We  focused 
on  compounds  involving  a  pentavalent  phosphorus 
atom  with  a  double  bond  with  oxygen,  the  common 
structure  of  the  species  of  interest. 

Melius  performed  ah  initio  calculations  with  the 
BAC-MP4  method  [20],  reporting  the  heat  of  for¬ 
mation  of  some  alkyl  phosphates  and  alkyl  phos- 
phonates  and  the  entropies  and  heat  capacities  of  six 
species  [21],  Taking  these  data  as  a  basis,  an  evalu¬ 
ation  of  group  additivity  contributions  was  used  to 
estimate  the  thermochemistry  of  a  large  number  of 
phosphorus  compounds.  First  six  groups  were  cal¬ 
culated.  Optimization  of  the  group  contributions  was 
performed,  using  as  a  first  guess  Benson's  values  [19] 
when  available  or,  by  analogy,  the  value  of  the  cor¬ 
responding  group  involving  (CO)  instead  of  (PO). 
Since  many  other  groups  are  needed,  an  evaluation 
of  the  properties  of  new  species  was  done  using  the 
UHF/PM3  Hamiltonian  and  the  MOPAC  program 
[22] .  Table  1  compares  heats  of  formation  calculated 
with  the  new  group  values,  BAC-MP4  method,  MO¬ 
PAC,  and  Bensons  group  values.  Despite  generally 
less  precise  PM3  evaluations,  the  difference  in  esti¬ 
mation  of  enthalpies  is  less  than  10  kcal  mol-1,  or 
6%,  below  the  estimated  error  given  by  Melius  [20], 
For  entropy  and  heat  capacity,  evaluations  done  by 
MOPAC  were  corrected  for  the  contribution  of  in¬ 
ternal  rotors  that  were  not  taken  into  account.  The 
difference  with  BAC-MP4  calculations  is  less  than  8 
cal  mol - 1,  or  8%,  for  entropy  and  appears  related  to 
groups  bonded  to  the  (PO)  group.  Correction  factors 
were  determined  and  applied  to  subsequent  calcu¬ 
lations.  Corrections  used  are,  for  each  substituent, 
—  1.5  cal  mol-1  for  OH,  —  1  cal  mol-1  for  CH3  or 
NH2,  —2  cal  mol-1  for  Q2H5,  —2.5  cal  mol-1  for 
OCH3,  and  +3.3  cal  mol-1  for  H,  F,  or  both. 

Nineteen  other  group  contributions  were  based 
on  evaluation  of  thermochemical  data  of  new  phos¬ 
phorus  compounds  by  MOPAC,  as  described,  in¬ 
volving  C,  II,  O,  F,  N,  S,  and  C  =  N  substituents. 


When  the  heat  of  formation  had  been  calculated  by 
Melius  [20],  the  BAC-MP4  value  was  preferred  for 
the  determination  of  the  enthalpy.  This  case  oc¬ 
curred  for  the  following  groups:  PO-(H)2(0),  PO- 
(H)(0)*  PO-(C)(F)*  PO-(C)(F)(0),  PO-(F)(0)* 
PO-(N)(0)*  and  C-(C)(H)2(PO).  Table  2  summa¬ 
rizes  all  proposed  group  values. 


Kinetic  Modeling 

Combustion  chemistry  of  organophosphorus  com¬ 
pounds  can  be  postulated  by  analogy  with  better- 
known  carbonated  species.  Radical  reactions  and 
molecular  eliminations  must  be  considered. 

Unimolecular  Decompositions 

Two  kinds  of  molecular  decomposition  seem  to 
play  a  role  in  the  degradation  of  alkyl  phosphate  or 
phosphonates.  A  six-center  concerted  reaction  leads 
to  elimination  of  an  alkyl  group,  producing  an  alkene 
and  an  OH  group  bonded  to  P.  The  reaction  occurs 
by  formation  of  a  cyclic  transition  state  where  an  H 
atom  of  the  alkyl  group  is  abstracted  by  the  O  atom 
of  the  (P  =  0)  [9,19].  The  four-center  reaction  in¬ 
volves  abstraction  of  H  atoms  from  an  OH  substit¬ 
uent  by  an  O  atom  that  leads  to  the  elimination  of 
alcohol  or  water  and  formation  of  a  PO,  group 

[9.19] . 

Zegers  and  Fisher  measured  rates  of  the  six-center 
elimination  for  DEMP  [9],  TEP  [10],  and  DIMP 
[11]  by  following  decay  of  the  reactant  in  pyrolysis 
between  700  K  and  900  K.  Rates  are  in  good  agree¬ 
ment  with  ester  decomposition  into  acid  and  alkene 
involving  the  same  kind  of  transition  state  [11].  How¬ 
ever,  measurements  of  ester  decomposition  show 
that  the  pre-exponential  factor  A  for  such  a  six-mem¬ 
ber  transition  state  is  between  3.1012  and  1013  s  1 

[11.19] ,  We  adopted  the  value  of  1013  s  - 1  for  DEMP 
decomposition,  considering  six  abstractable  H  atoms 
instead  of  three.  Because  of  the  larger  H  atom  num¬ 
ber,  A  should  be  about  two  times  higher  for  the  re¬ 
action  of  TEP  and  DIMP.  Activation  energies  were 
corrected  to  fit  these  A  factors  to  experimentally 
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TABLE  2 

Group  contributions  for  organophosphorus  compounds  at  298  K 


Group 

AH®  (kcal  mol  x) 

S°  (cal  K  1  mol  b) 

Cb  (cal  K  1  mol  b) 

P0-(C)(H)(0) 

-74.3 

23.1 

9.5 

P0-(C)(0)2 

-92.7 

0.3 

9.3 

PO-(0)3 

- 106.4 

-0.2 

9.1 

PO-(H)2(0) 

-69.1 

39.7 

9.5 

P0-(H)(0)2 

-91.1 

19.2 

9.3 

PO-(C)(F)(H) 

- 138.0 

49.3 

12.1 

PO-(C)(F)2 

-215.1 

55.3 

16.5 

P0-(C)(F)(0) 

-161.6 

24.7 

11.4 

P0-(F)(0)2 

- 173.1 

20.6 

12.1 

PO-(F)(H)(0) 

- 157.5 

47.7 

12.6 

PO-(CN)(C)(0) 

-46.3 

30.9 

15.4 

P0-(C)(N)(0) 

-  70.3 

4.8 

9.9 

P0-(N)(0)2 

-94.5 

-0.5 

9.6 

P0-(CN)(N)(0) 

-42.0 

33.7 

17.1 

P0-(C)(0)(S) 

-63.9 

1.6 

6.9 

C-(H)3(PO)a 

-9.95 

28.1 

6.19 

C-(C)(H)2(PO) 

-5.2 

3.0 

3.9 

0-(H)(PO) 

-54.9 

27.9 

4.7 

0-(C)(PO) 

-36.8 

8.2 

2.8 

0-(PO)2 

-.56.4 

1.8 

6.2 

N-(H)2(PO)b 

-14.9 

24.69 

4.07 

N-(C)(H)(PO) 

-3.2 

4.4 

3.1 

N-(C)2(PO) 

8.6 

-17.3 

4.5 

S-(H)(PO)° 

-1.4 

31.2 

7.63 

S-(C)(PO) 

6.4 

11.3 

5.4 

aC-(H)3(PO)  =  C-(H)3(CO)  assigned  for  C°. 
bN-(H)2(PO)  =  N-(H)2(CO)  assigned. 
°S-(H)(PO)  =  S-(H)(CO)  assigned. 


measured  rate  constants  at  800  K.  In  the  case  of 
TEP,  A  is  higher  to  keep  the  same  activation  energy 
as  for  DEMP.  The  rates  are  given  in  Table  3. 

No  experimental  measurements  have  been  re¬ 
ported  for  the  four-center  decomposition.  To  esti¬ 
mate  these  rates,  experimental  results  of  Zegers  and 
Fisher  were  modeled.  Based  on  the  analysis  of  re¬ 
action  products,  alkene  and  alcohol,  we  assumed, 
like  these  authors,  that  no  radical  reaction  occurs 
with  the  parent  compound,  but  only  successive  six- 
center  and  four-center  eliminations.  The  postulated 
mechanism  is  given  in  Table  3.  Rates  of  the  six-cen¬ 
ter  decomposition  of  the  intermediates  were  the 
same  as  the  reactant  decomposition  after  a  correc¬ 
tion  of  the  A  factor  for  taking  into  account  the  H 
atom  number.  Rates  of  the  four-center  eliminations 
were  then  adjusted  to  reproduce  experimentally  ob¬ 
served  formation  of  alcohol  and  alkene. 

Experimental  results  [9-11]  were  obtained  in  an 
atmospheric  pressure  flow  reactor,  with  the  reactant 
diluted  in  nitrogen.  Measured  products  were  ethyl¬ 
ene  and  ethanol  in  the  reaction  of  DEMP  and  TEP, 
and  propene  and  1-propanol  in  the  reaction  of 


DIMP.  Because  heterogeneous  effects  change  the 
ratio  between  the  products  [9-11],  the  initial  load 
used  for  the  calculations  was  the  amount  extrapo¬ 
lated  to  t  —  0  at  the  lower  temperature  studied.  This 
assumption  supposes  that  heterogeneous  reactions 
are  independent  of  temperature  and  play  a  role  in 
the  entrance  section  of  the  flow  reactor  but  may  be 
neglected  downstream.  For  DIMP,  the  load  used 
was  exactly  the  initial  one,  but  the  start  of  the  reac¬ 
tion  was  shifted  by  10  ms. 

Pre-exponential  factors  for  the  four-center  elimi¬ 
nations  were  chosen  to  be  in  the  same  order  as  cor¬ 
responding  reactions  reported  with  carbonated  com¬ 
pounds  [19],  respectively,  2.5  X  1013,  3.5  X  1013, 
and  5  X  1013  s_1  for  formation  of  ethanol,  i-pro- 
panol,  and  water.  With  such  values,  the  activation 
energy  should  be  near  the  enthalpy  of  reaction  to  be 
able  to  reproduce  the  experimental  results,  a  situa¬ 
tion  veiy  different  from  other  known  reactions  in¬ 
volving  the  same  kind  of  transition  state  [18].  The 
activation  energy  of  the  reverse  reaction,  a  molecular 
addition  of  water  or  alcohol  on  a  POa  function,  is 
then  very  low  and  should  be  below  2  kcal  mol  ~~ 1 . 
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TABLE  3 

Pyrolysis  mechanism  of  DEMP,  DIMP,  and  TEP 


Reaction 

A  (s' 

-1) 

b 

E  (cal  mol  4) 

DEMP  PO(OH)(OEt)(CH3)  +  C2H4 

1 

X 

1013 

0 

45,300“ 

PO(OH)(OEt)(CH3)  <b4  PO(OH)2(CH3)  +  C2H4 

5 

X 

1012 

0 

45,300 

PO(OH)(OEt)(CH3)  <H>  CH3P02  +  C2H5OH 

2.5 

X 

1013 

0 

44,000 

PO(OH)2(CH3)  <h.  ch3po2  +  h2o 

5 

X 

1013 

0 

39,000 

DIMP  «  PO(OH)(OiPr)(CH3)  +  C3H6 

2 

X 

1013 

0 

41,200“ 

PO(OH)(OiPr)(CH3)  o  PO(OH)2(CH3)  +  C3H6 

1 

X 

1013 

0 

41,200 

PO(OH)(OiPr)(CH3)  >H>  CH3P()2  +  C3H7OH 

3.5 

X 

1013 

0 

43,000 

TEP  e->  PO(OH)(OEt)2  +  C2H4 

2.8 

X 

1013 

0 

45,300“ 

PO(OH)(OEt)2  PO(OH)2(OEt)  +  C2H4 

1.9 

X 

1013 

0 

45,300 

PO(OH)(OEt)2  <H>  C2H50P02  +  C2HsOH 

2.5 

X 

1013 

0 

44,000 

PO(OH)2(OEt)  e->  PO(OH)3  +  C2H4 

9.6 

X 

1012 

0 

45,300 

PO(OH)2(OEt)  HOPOo  +  C2H5OH 

2.5 

X 

1013 

0 

44,000 

c2h5opo2  hopo2  +  C2H4 

1.9 

X 

1013 

0 

45,300 

PO(OH)2(OEt)  e->  C2H50P02  +  H20 

5 

X 

1013 

0 

45,000 

PO(OH)3  o  hopo2  +  h2o 

1.5 

X 

1014 

0 

46,000 

“Calculated  from  Refs.  [9-11],  see  text. 


Residence  time  (s) 


Fig.  2.  Profiles  of  species  in  the  pyrolysis  of  DIMP  at 
799  K  in  flow  reactor  under  1  atm  [11],  Curves  are  model 
predictions. 


Figure  1  compares  the  calculated  ratios  of  alcohol 
and  alkene  to  the  experimental  ones  for  the  three 
reactants.  For  the  highest  conversions,  obtained  for 
long  residence  time,  the  alcohol  production  drops 
off  in  agreement  with  experimental  results.  This  is 
due  to  a  shift  in  the  equilibrium  of  the  four-center 
elimination  reactions,  which  are  in  partial  equilib¬ 
rium  at  high  conversion.  This  shift  favors  the  reverse 
reaction  of  addition  of  alcohol. 


Fig.  1.  Ratio  of  alcohol  and  alkene  amount  formed  in  CH3P02  +  C2HsOH  PO(OH)(OEt)(CH3) 
the  pyrolysis  of  DEMP  [9],  TEP  [10],  and  DIMP  [11]  in 

flow  reactor  under  1  atm.  Curves  are  model  predictions.  which  produces  the  olefin, 
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Fig.  3.  Experimental  and  computed  profiles  of  the  tem¬ 
perature  in  the  stoichiometric  premixed  H2/02/Ar  flame 
[13]  (mole  fraction:  0.26/0.13/0.61),  doped  with  0.2%  vol. 
of  DMMP.  P  =  80  Torr. 


Fig.  4.  Experimental  and  computed  profiles  of  H2,  02, 
and  H20  in  the  stoichiometric  premixed  H2/02/Ar  flame 
[15]  (mole  fraction:  0.26/0.13/0.61)  with  and  without  the 
addition  of  0.2%  vol.  of  TMP  P  =  47  Torr.  Curves  are 
model  predictions.  Filled  symbols  and  curves  refer  to  the 
doped  flame,  open  symbols  and  dashed  curves  refer  to  the 
undoped  flame. 

PO(OH)(OEt)(CH3)  PO(OH)2(CH3)  +  C2H4 

Fig.  2  shows  reactants  and  carbonated  products  pro¬ 
files  for  the  pyrolysis  of  DIMP.  The  most  important 
phosphorus  product  is  then  CPI3P02,  which  is  in 
equilibrium  with  water  and  i-propanol. 

Combustion  Model 

A  mechanism  has  been  built  to  reproduce  com¬ 
bustion  of  DMMP  and  TMP.  It  includes  reaction 


types  usually  considered  in  the  combustion  of  hy¬ 
drocarbons  and  the  possible  molecular  eliminations. 
Major  molecular  intermediates  taken  into  account 
are  PO(OH)(OCH3)2,  PO(OH)(CH3)(OCH3), 
PO(OH)2(OCH3),  PO(OH)2(CH3),  PO(OH)3, 
CH3PO,  CH3OPO,  CH3P02,  and  CH30P02.  These 
compounds  react  by  molecular  elimination,  bond- 
scission,  II  abstraction  by  02  or  a  free  radical,  and 
addition  of  H  or  OH  to  the  P  =  O  double  bond.  The 
radicals  produced  then  react  by  combination  with 
small  free  radicals  or  by  decomposition  involving  a 
bond  breaking  via  /?-scission.  The  ultimate  products 
are  HOPO,  HOP02,  I  IPO,  HPOz,  PO,  P02,  and 
P03.  For  these  species,  all  unimolecular  and  bimo- 
lecular  possible  reactions  are  written  by  considering 
also  reactions  with  H,  O,  OH,  HOa,  H20,  and  CH3. 
The  complete  mechanism  consists  of  41  different 
phosphorus  species  involving  202  reactions.  The  ki¬ 
netic  data  were  estimated  by  analogy  with  reactions 
of  carbon-  or  nitrogen-containing  species  and  by 
thermochemical  considerations.  The  molecular 
elimination  rates  were  evaluated  by  analogy  with 
ones  determined  in  the  pyrolysis  study.  Rates  of  H 
abstraction  were  calculated  with  the  Evans-Polanyi 
correlation  proposed  by  Dean  and  Bozzelli  [23],  The 
addition  of  H  and  OH  to  P  =  O  were  estimated,  re¬ 
spectively,  at  5.0  X  1012  exp(  —  5000 !RT)  cm3  mol-1 
s-1  and  1.4  X  1012  exp(  -  2000/RT)  cm3  mol“ 1  s “ 1 
by  analogy  with  the  addition  to  a  C  =  O  bond. 

The  mechanism  used  for  modeling  the  combus¬ 
tion  of  oxygenated  and  carbonated  species  was  de¬ 
rived  from  a  mechanism  of  combustion  of  dimeth- 
ylether  [24] .  We  retained  species  containing  two  or 
fewer  C  atoms  and  their  associated  reactions,  giving 
a  submechanism  of  51  species  and  263  reactions. 

This  mechanism  has  been  used  to  model  the  com¬ 
bustion  of  DMMP-  and  TMP-doped  H2/02/Ar  sta¬ 
bilized  premixed  stoichiometric  flames  in  the  exper¬ 
imental  conditions  of  Korobeinichev  [13,15],  The 
calculations  were  performed  with  the  PREMIX  [25] 
code.  Fig.  3  compares  computed  and  experimental 
temperature  profiles  of  the  flame  with  and  without 
the  addition  of  0.2%  DMMP  [13].  The  increased 
heat  release  when  DMMP  is  added  is  well  repro¬ 
duced  by  the  model.  Fig.  4  compares  the  profiles  of 
reactants  and  water  obtained  with  and  without  the 
addition  of  TMP  at  the  experimental  conditions  of 
Ref.  [15],  For  both  DMMP  and  TMP,  the  promoting 
effect  of  the  phosphorus  additive  is  well  reproduced 
by  the  model.  The  reaction  mechanism  is  available 
from  the  authors  or  on  the  Web  at  http://www- 
cms.llnl.gov/combustion/combustion_home.html. 

The  additives  are  quickly  consumed  in  the  first 
stage  of  the  flame.  Intermediate  alkyl  phosphates 
and  phosphonates  are  produced  close  to  the  burner 
by  addition  of  I  I  to  the  P  =  O  bond,  followed  by  de¬ 
composition. 
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Fig.  5.  Computed  profiles  of  DMMP  and  intermediate 
phosphorus  species  in  the  stoichiometric  premixed  H2/02/ 
Ar  flame  [13]  (mole  fraction:  0.26/0.13/0.61),  doped  with 
0.2%  vol.  of  DMMP,  P  =  80  Torr. 


Fig.  6.  Experimental  and  computed  profiles  of  HOPO, 
HOP02,  and  orthophosphoric  acid  in  the  stoichiometric 
premixed  H2/02/Ar  flame  [15]  (mole  fraction:  0.26/0.13/ 
0.61)  doped  with  0.2%  vol.  of  TMP.  P  =  47  Torr.  Curves 
are  model  predictions. 

PO(CH3)(OCH3)2  +  II 
->  PO(OH)(CH3)(OCH3)  +  ch3 

Fig.  5  shows  computed  profiles  of  DMMP  and  phos¬ 
phorus  intermediates  at  the  experimental  conditions 
of  Ref.  [13].  These  compounds  react  mainly  by  H 
atom  addition  or  a  four-center  elimination 

PO(OH)(CH3)(OCH3)  +  II 
-»  PO(OH)2(CH3)  +  ch3 
PO(OH)(CH3)(OCH3)  ch3po2  +  ch3oh 


Fig.  7.  Experimental  and  computed  profiles  of  PO  and 
P02  in  the  stoichiometric  premixed  H2/02/Ar  flame  [15] 
(mole  fraction:  0.26/0.13/0.61)  doped  with  0.2%  vol.  of 
TMP.  P  =  47  Torr.  Curves  are  model  predictions. 

When  temperature  increases,  the  H-abstraction 
channel  becomes  predominant  and  leads  to  succes¬ 
sive  decompositions: 

PO(CH3)(OCH3)2  +  H 
-4  PO(CH3)(OCH3)(OCH2)  +  H2 

PO(CH3)(OCH3)(OCH2) 

-4  PO(CH3)(OCH3)  +  ch2o 

PO(CH3)(OCH3)  -4  CH3OPO  +  ch3 
By  addition  of  H  radical,  species  as  CH3P02  or 
CH3OPO  produce  small  oxides.  These  are  ortho¬ 
phosphoric  acid,  HOPO,  HOP02,  and  the  radicals 
PO  and  P02.  Figs.  6  and  7  compare  experimental 
and  calculated  profiles  of  these  species  in  the  flame 
doped  with  TMP  [15],  Near  the  burner,  orthophos¬ 
phoric  acid  is  most  stable,  decomposing  into  H0PO2 
and  water  when  the  temperature  increases.  HOPO 
reacts  mainly  by  H  abstraction,  giving  P02,  while  H 
adds  to  HOP02  to  give  PO(OH)2  radical,  which  then 
combines  with  other  radicals  at  lower  temperatures 
or  decomposes  into  HOPO  and  OH  at  higher  tem¬ 
peratures. 

The  promoting  effect  of  the  phosphorus  additive 
comes  from  enhanced  heat  release  closer  to  the 
burner  from  exothermic  recombination  reactions  in¬ 
duced  by  the  phosphorus  species.  In  other  condi¬ 
tions,  such  as  methane  flames  [1-4],  the  drop  of  rad¬ 
ical  concentration  in  the  flame  caused  by 
recombination  leads  to  an  inhibiting  effect:  in  the 
case  of  the  H2/02  flame,  where  II  radicals  are  veiy 
abundant  and  react  mainly  with  oxygen,  the  higher 
temperature  allows  a  strong  acceleration  of  the  com¬ 
bustion  by  the  way  of  the  very  sensitive  branching 
reaction 
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h  +  o2^oh  +  o 

with  a  high  activation  energy  of  16.5  kcal  mol-1. 

The  sensitivity  analysis  shows  that  in  the  very  be¬ 
ginning  of  the  H2/02-doped  flame,  a  termination  cy¬ 
cle  appears  with  H  abstractions  on  the  DMMP  or 
TMP  followed  by  combination  of  the  produced  rad¬ 
ical  or  the  radicals  produced  after  a  decomposition. 
In  the  case  of  DMMP, 

PO(CH3)(OCH3)(OCH2)  +  H^PO(CH3)(OCH3)2 

Later,  the  major  termination  effect  conies  from  the 
cycle  already  proposed  by  Twarowski  [18]  that  in¬ 
volves  P02: 

P02  +  H  — >  HOPO 

HOPO  +  X  ^  P02  +  HX 

where  X  =  H,  O,  or  OH.  Despite  the  loss  of  H 
atoms,  which  react  otherwise  with  02  and  produce 
chain  branching,  these  termination  cycles  have  a 
promoting  effect  because  of  the  high  exothermicity 
of  H2  or  water  formation.  The  present  modeling  re¬ 
sults  indicate  that  formation  of  the  more  stable 
H0P02,  which  reduces  the  termination  cycles,  has 
an  inhibiting  effect  by  reducing  the  amount  of  P02 
radicals: 

P02  +  OH  ^  H0P02 

In  contrast,  the  consumption  of  this  species  to  give 
HOPO  increases  the  global  reaction  rate: 

H0P02  +  H  -A  PO(OH)2 

PO(OH)2  -a  HOPO  +  OH 

Conclusion 

An  evaluation  of  thermodynamical  properties  of 
organophosphorus  compounds  was  performed  by 
proposing  some  incremental  group  values  for  en¬ 
thalpy,  entropy,  and  heat  capacity.  The  25  groups  in¬ 
volving  C,  H,  N,  O,  F.  S,  and  CN  bonded  to  PO 
allow  prediction  of  properties  of  many  species,  for 
which  no  data  were  currently  available,  with  reason¬ 
able  accuracy. 

Modeling  of  pyrolysis  of  DEMP,  TEP,  and  DIMP 
permitted  evaluation  of  the  rate  constants  of  some 
molecular  eliminations  whose  importance  was  pre¬ 
viously  not  easy  to  estimate  under  combustion  con¬ 
ditions.  A  large  combustion  mechanism  developed 
for  DMMP  and  TMP  was  validated  against  experi¬ 
mental  results  obtained  in  doped  H2/02  flames.  This 
new  mechanism,  taking  into  account  all  possible  in¬ 
termediates,  is  more  complete  and  comprehensive 
than  previous  ones.  The  model  allows  the  prediction 
of  decomposition  routes  and  the  products  formed, 
which  are  interesting  in  incineration  of  chemical 


waste.  An  explanation  of  the  promoting  effect  on  the 
studied  flames  is  proposed  and  evaluated  for  each  of 
the  different  phosphorus  oxides. 
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COMMENTS 


Richard  G.  Gann,  NIST,  USA.  I  am  glad  to  see  this 
mechanistic  approach.  Our  work  in  the  United  States 
next  generation  fire  suppression  technology  program  in¬ 
dicates  that  chemically,  phosphorus-containing  com¬ 
pounds  are  effective  in  proportion  to  the  number  of  P 
atoms  added  to  the  flame.  This  work  is  consistent  with 


the  experimental  finding.  The  exception  may  be  the  cy¬ 
clic  phosphonitriles,  in  which  the  phosphorus  is  very 
tightly  bound.  The  Army  Research  Laboratory  is  looking 
into  this  at  present. 

Author’s  Reply.  The  comment  requires  no  reply. 


